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ABSTRACT 

The  e'fect  of  particles,  such  as  dust  in  air,  on  aerodynamic 
drag  of  circular  cylinders  was  calculated  for  compressible  flow  at 
critical  Mach  number  and  for  incompressible  flow.  The  effect  of 
compressibility  was  found  negligible  for  particles  larger  than  about  10 
pm,  for  which  the  air  can  be  considered  a  continuum.  Drag  coefficient 
and  collection  efficiency  are  provided  for  a  wide  range  of  inertia 
parameters  and  Reynolds  numbers  for  both  compressible  and 
incompressible  flow. 
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NOTATION 

d  particle  diameter,  cm 

Em  impaction  efficiency  of  particles  on  a  target 

F  force  exerted  on  the  cylinder  by  the  particles  alone  travelling 

in  air,  dynes 

L  cylinder  radius,  cm 

p  drag  pressure  on  cylinder  due  to  stream  of  particles  alone, 

dynes  cm"2 

r  distance  from  any  point  to  the  origin,  cm 

t  time,  seconds 

U  free-stream  velocity,  cm  s"1 

u  local  fluid  velocity,  cm  s"1 

ux  parallel  component  of  fluid  velocity,  cm  s"1 
uy  transverse  component  of  fluid  velocity,  cm  s"1 
v  local  particle  velocity,  cm  s'1 

vx  component  of  particle  velocity  parallel  to  free  field  flow 
direction  immediately  before  impact,  cm  s"1 

vx'  component  of  particle  velocity  parallel  to  flow  direction 
immediately  after  impact  and  reflection,  cm  s"1 
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far  upstream  particle  velocity,  cm  s'* 


component  of  particle  velocity  perpendicular  to  free  field  flow 
direction  immediately  before  impact,  cm  s'1 


change  in  the  component  of  particle  velocity  parallel  to  the 
free  stream  flow  direction  caused  by  reflection  from  the  sphere 
or  cyl inder ,  cm  s' 1 


co-ordinate  (origin  at  centre  of  sphere  or  cylinder)  of 
particle  position  parallel  to  free  stream  flow  direction,  cm 


transverse  co-ordinate  of  particle  position,  cm 


angle  of  incidence  of  a  particle  against  the  sphere  or 
cylinder,  radians 


angle  between  particle  direction  and  free  field  flow  direction 
just  before  the  particle  hits  the  sphere  or  cylinder,  radians 

total  angle  change  of  a  particle  from  its  far  upstream 
direction  to  its  direction  after  reflection  from  the  cylinder 
or  sphere,  radians 


polar  angle  between  negative  x  axis  and  radius  vector  to 
particle  impact  position,  radians 


mean  free  path  of  molecules  pm 


local  absolute  viscosity  of  fluid,  poise 


free  stream  absolute  viscosity  of  fluid,  poise 
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local  fluid  density,  g  cm-3 


fluid  density  at  free  stream  Mach  number  g  cm'3 


fluid  density  of  still  air  g  cm'3 

uniform  particle  mass  per  unit  volume  of  air  far  upstream,  g 
cm'3 


particle  density  g  cm'3 

angle  between  positive  x  axis  and  position  in  flow  field 
radians 

The  following  were  dimensionless: 

Cunningham  Correction  Factor 


drag  coefficient  due  to  particle  impact  with  cylinder 


change  in  component  of  particle  velocity  due  to  reflection  from 
sphere  or  cylinder  as  a  function  of  off-axis  distance  far 
upstream 


particle  inertia  parameter 


Mach  number 


spherical  particle  Reynolds  number  in  flow  influenced  by 
presence  of  cylinder 


spherical  particle  Reynolds  number  in  free  stream 
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r 


1 


u 


v 


V 

X 


r,  distance  from  any  point 
L 


time 


tu 


local  fluid  velocity  uv 

U  x 


local  particle  velocity  jj, 


dx 

dl 


,  parallel  component  of 


to  the  origin 


and  Uy  are  x  and  y  components 

vY  and  v..  are  x  and  y  components 

*  y 


particle  velocity 


v 


x 


y 


y 

y 


Vd 

*i 

ra 

0 


transverse  component 
dT 


parallel  co-ordinate  ^ 


of  particle  velocity 


y 

transverse  co-ordinate  - 

L 

off-axis  distance  of  a  particle  at  point  of  impact  with  the 
cyl inder 


far  upstream  transverse  co-ordinate  of  the  envelope  of 
particles  which  eventually  hit  the  cylinder 

co-ordinate  of  particle  used  in  discretization  for  integration; 
y^  increases  with  i  from  yi  =  o  to  yn  =  yd 


specific  heat  ratio  for  air 


dimensionless  group  independent  of  drop  size  formed  by  combin¬ 
ing  Re„  and  K 
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INTRODUCTION 


Extensive  studies  have  been  made  of  the  drag  exerted  on  trans¬ 
verse  cylinders  in  low  speed  steady  flow  of  moving  air  for  which  com¬ 
pressibility  effects  are  negligible  [1],  These  results  are  shown  for 
Reynolds  numbers  up  to  10s,  which  barely  extends  into  the  supercritical 
Reynolds  number  region.  Further  work  with  large  cylinders  in  a  wind 
tunnel,  with  velocities  low  enough  so  that  compressibility  effects  are 
still  negligible,  have  provided  information  about  aerodynamic  drag 
reaching  further  into  the  supercritical  region  with  Reynolds  numbers 
from  10s  to  107  [2],  Experiments  in  which  the  flow  velocities  are  high 
enough  to  provide  an  effect  on  drag  of  cylinders  due  to  compressibility 
have  also  been  performed  [3].  Similar  information  has  been  provided 
for  Spheres  [1],  Some  information  is  also  available  for  unsteady  in¬ 
compressible  flow  which  was  obtained  from  field  experiments  with  large 
blast  waves  [4,  5].  The  reports  cited  generally  show  that  drag  coeffi¬ 
cient  increases  with  Mach  number  for  Mach  numbers  between  0.3  and  1. 
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The  blast  waves  from  500  ton  TNT  bursts  on  the  DRES  test  range 
were  associated  with  visible  quantities  of  dust  picked  up  from  the  dry 
prairie  ground.  This  was  observed  by  the  author  during  Operation 
Prairie  Flat  in  1968  and  in  Event  Dial  Pack  in  1970.  Therefore,  to 
fully  understand  the  forces  on  targets  in  such  a  blast  wave,  a  know¬ 
ledge  of  the  effect  of  the  dust  is  necessary.  Attempts  with  laboratory 
experiments  hav*'  been  performed  to  determine  the  effect  of  introducing 
dust  into  the  air  stream  [6,  7]  and  to  calculate  the  effect  using  a 
mathematical  model  [8].  Only  incompressible  flow  was  considered  in 
these  experimental  and  theoretica"  studies. 

The  purpose  of  this  report  is  to  describe  numerical  solutions  to 
the  equations  of  motion  of  particles  moving  toward  and  impacting  with  a 
transverse  circular  cylinder  in  compressible  air  flow.  The  solutions 
are  subsequently  applied  to  calculate  drag  coefficients  due  to  parti¬ 
cles  alone  and  tha  significance  of  compressibility  is  described. 

DEFINITION  OF  THE  PROBLEM 

A  cloud  of  spherical  particles  uniformly  distributed  in  steady 
air  flow  is  assumed  to  travel  toward  a  circular  cylinder  with  its  axis 
perpendicular  to  the  flow  direction.  Far  upstream,  the  particles  are 
assumed  to  travel  along  with  the  velocity  and  direction  of  the  air,  but 
as  they  approach  the  cylinder  they  will  not  necessarily  follow  stream¬ 
lines  because  the  air  velocity  changes  in  both  speed  and  direction  in 
the  vicinity  of  the  cylinder.  Far  upstream,  the  particles  which  travel 
toward  the  cylinder  and  eventually  contact  it  will  be  contained  within 
a  width  less  than  the  cylinder  diameter.  This  distance  or  limiting 
envelope  divided  by  the  cylinder  diameter  is  called  the  collection 
efficiency.  The  velocity  of  a  particle  immediately  before  contact  and 
also  the  point  of  contact  with  the  cylinder  depend  upon  the  position  of 
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the  particle  within  the  limiting  envelope  far  upstream  from  which  it 
travels  to  the  target.  All  particles  have  their  momentum  changed  when 
they  collide  with  the  cylinder  and,  in  the  case  of  solid  particles,  are 
deflected  in  a  new  direction.  The  total  change  in  momentum  of  all 
particles  which  hit  the  cylinder  is  balanced  by  the  force  they  exert  on 
it.  The  basic  problem  is,  therefore,  to  calculate  the  collection 
efficiency  and  then  to  calculate  the  force  of  particles  on  the 
cylinder,  from  which  the  drag  coefficient  due  to  the  particles  alone 
can  be  obtained. 

EQUATIONS  OF  MOTION 


The  motion  of 

an  individual  spherical  particle  has 

been  shown 

[9]  to 

be  determined 

by  the  following  ordinary  differential 

equations : 

dv 

.  C0Re(i}y  -  V 

24K 

dT 

(1) 

d\ 

CDRe(“x  -  V 

dT 

24K 

(2) 

where 

Re 

=  Re°  [(%  -  vy)2  +  (ux  -  vx)2]? 

(3) 

K 

0(j2\j 

=  j-  particle  inertia  parameter 

(4) 

Re0 

_  UdRo  free  stream  Reynolds  number 

Mo 

(5) 

The  symbols  are  defined  in  the  Notation  Section  near  the  front  of  this 
report  and  the  basic  geometry  of  the  flow  system  is  illustrated  in  Fig. 
1. 
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The  above  equations  of  motion  are  for  incompressible  flow.  In 
compressible  flow  the  density  changes  as  flow  velocity  changes  near 
the  cylinder.  Then  equation  (3)  must  be  replaced  by  the  following 
equation: 


Re  =  Mkr)[<0y- V2  ^x-^x)2]* 

The  density  ratio  p/p0  is  calculated  from  the  energy  equation  for 
isentropic  flow  [10]  which  is  written  as  follows: 


i  ♦  yi  Mo2 
2 


1 


V1 


(7) 


where  ps  is  for  M0  =  0  and  p0  is  for  free  stream  Mach  number,  M0 . 


A1  so 


£s 

P 


1  +  Ta-1  M2 


1 

\  T,-l 


where  p  is  for  the  local  Mach  number  M,  near  the  cylinder  which 
given  by  the  following  equation: 


(8) 

is 


M 


Mn^U  2  +  U  2 


(9) 


Then  by  combining  equations  (7),  (8)  and  (9),  the  ratio  of  densities  in 
equation  (6)  is  obtained  in  terms  of  the  free  stream  Mach  number  as 
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f o 1  lows : 


P 

Po 


1  +  i  (V1)  M°2 

l  +  i  (ra-i)  M02  (Uy2  +  Uy2) 


(10) 


Several  assumptions  were  inherent  in  the  development  of  equation 
(1)  and  equation  (2)  for  calculating  impaction  efficiency  and  force  due 
to  a  stream  of  particles,  including: 


a.  uniform  particle  distribution; 

b.  no  gravitational  or  electrostatic  forces  of  consequence; 

c.  monodisperse  spherical  particles  with  diameter  negligible 
in  relation  to  the  target  sphere  or  cylinder  diameter;  and 

d.  free  stream  flow  that  was  steady  and  irrotational ; 

e.  the  ratio  of  the  local  viscosity  p  to  the  free  stream 
viscosity  p0,  is  unity.  This  assumption  is  valid,  because 
viscosity  is  not  appreciably  affected  by  the  order  of 
magnitude  of  the  changes  in  temperature  found  along  a 
trajectory. 


The  drag  coefficient  is  a  function  of  Reynolds  number  and  was 
available  in  the  form  of  definitive  empirical  equations  [11].  These 
equations  are  stated  as  follows: 
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Re 


'-'rrc 

-  2.3363  x  10~4(CgRe2)2  +  2.0154  x  10‘6(CDRe2)3 
6.9105  x  10-9(C0Re2)" 


for  Re  <  4  or  CpRe2  <  140 


(11) 


logloRe  =  -1.29536  +  9.86  x  10"1  (logloCDRe2)  -  4.6677  x 

10’2(1ogloC0Re2)2  +  1.1235  x  10'3  (log10CDRe2)3  (12) 

for  3  <  Re  <  104  or  100  <  CQRe2  <  4.5  x  107 

Equation  (11)  approaches  Stokes  flow  as  the  Reynolds  number  approaches 
zero.  These  equations  assume  that  the  fluid  is  a  continuum.  For 
Knudsen  number,  which  is  the  ratio  of  the  molecular  mean  free  path  in 
the  fluid  to  the  diameter  of  the  particle,  less  than  a  value  of  the 
order  of  10'3,  the  particle  is  large  by  comparison  with  the  scale  of 
molecular  processes  in  the  continuous  phase  and  the  fluid  can  be 
treated  as  a  continuum.  The  mean  free  path  for  air,  at  293  K  tempera¬ 
ture,  and  pressure  of  100  kP,  is  0.0653  pm.  Thus,  noncontinuous  ef¬ 
fects  are  significant  for  particles  smaller  than  about  10  pm  [12].  The 
actual  significance  on  collection  efficiency  and  drag  coefficient  was 
estimated  by  dividing  the  drag  coefficient  given  by  equations  (11)  and 
(12)  by  the  Cunningham  correction  factor  and  calculating  the  effects 
for  some  of  the  smaller  particle  sizes  considered.  The  factor  is  ap¬ 
plied  to  account  for  the  fact  that  the  fluid  resistance  is  smaller  than 
that  given  by  Stokes  law  when  noncontinuum  effects  are  present.  The 
correction  factor  is  written  as  follows  [13]: 

2\  r  ,-l.ldx-, 

Cr  =  1  +  —  1 1 .257  +  0.4  exp  ( - )] 

c  d  1  K  V  2X  n  (13) 

AIR  FLOW  FIELD 


The  compressible-air  flow  field  around  a  right  circular  cylinder 
of  infinite  extent  for  flow  perpendicular  to  the  cylinder  axis  is 
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described  by  the  following  equations  [14]: 


M2 


-  r"2)  cos2u  +  (1 

+  F*2) 

sin2u>]  + 

f _  1 3_p - 2  +  3--i,  _ 
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I?-  - 

-I?-6) 
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12  > 

v4 
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16 
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+  -i?-4 

*  -  If- 
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{  48 
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24 
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) 
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16 
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[-r'2  sin  2  toj  + 
M2 
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4 
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-  — r~6 
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\ 

-  — r'* 
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80 
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- 5f-2  _  _5f-<*  +  If-*  sin  5  u  Cos  u 

16  16  4 


sin  3  to  cos  to  + 
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The  first  bracketed  term  describes  the  incompressible  flow  around  the 
cylinder.  The  term  with  the  M2  factor  and  the  term  with  the  M*1  factor 
are,  respectively,  the  first  and  second  order  corrections  for  the  ef¬ 
fect  of  compressibility.  In  equations  (14a,  b)  u>  is  measured  counter¬ 
clockwise  from  the  x  -  axis  of  the  cylinder  as  is  shown  in  Figure  1. 
Equations  (14a,  b)  a^e  not  the  same  as  the  equations  stated  in 
reference  [14],  but  close  inspection  reveals  that  the  difference  is  due 
only  to  an  er^or  in  writing  the  equations  in  the  reference.  Equations 
(14a,  b)  were  derived  from  a  velocity  potential  function  [15], 

DRAG  FORCE  DUE  TO  PARTICLES  ALONE 

The  force  F^  of  a  stream  of  particles  hitting  a  cylinder  is 
calculated  from  the  change  in  momentum.  Assuming  perfectly  elastic 
spherical  particles,  that  the  particles  do  not  interact  with  each 
other,  and  that  only  one  impaction  of  each  particle  occurs,  or  at  least 
is  significant,  the  drag  force  is  given  by  the  following  equation: 

F  ^d  _ 

FTFT=  2V,/4Y*  (15) 

Avx  is  the  change  of  the  particle  velocity  component  in  the  direction 
of  free  stream  flow  due  to  reflection  from  the  cylinder.  It  is  given 
by: 


Av  =  V  -  V  COS  1 
X  X 


(16) 


where 


>v2 

x 


+  V 


(17) 
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and  1  is  the  angle  the  reflected  particle  makes  with  the  free  stream 
flow  direction  as  shown  in  Figure  2. 


r 


n 


2  tan'1 


- 1 
y 

-  tan'1 

(V 

JL 

'  y/  _  1 2  1 

i-y 


The  drag  coefficient  is  defined  as  follows: 


(18) 


(19) 


where  p  is  the  average  drag  pressure  [1]. 

Then,  using  equations  (15),  (16),  and  (19),  the  drag  coefficient  is 
given  by  the  following  equation: 

CD  =  2\o  P  '  *  cos  r>d*  (20) 

0 

A  more  detailed  derivation  of  equation  20  is  shown  elsewhere  [8]. 

Equation  20  was  solved  numerically  by  dividing  the  upstream 
cross  sectional  area  into  strips  of  equal  width  Ay,  and  calculating  the 
contribution  from  each  strip.  These  contributions  were  then  summed  as 
fol lows : 


Let  vy  -  v  cos  1  -  f(y) 


(21) 
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Then  the  integral  in  equation  (20)  is  written: 


yd  _  _  *  f (y i )  +  f(yn)  _  i=n-l  _  _ 

/  f(y)dy  =  - -  Ay  +  l  f(y,)Ay  (22) 

o  2  i=2  1 

where  y.  increases  with  i  from  yi  =  o  to  y  =  y  ,. 
i  n  o 

To  solve  equation  (22)  to  obtain  the  drag  coefficient  given  by  equation 

(20)  it  is  necessary  to  know  the  values  of  v  ,  v,  and  Y  in  the  range 

0  -<  y  ^  Jd. 

The  first  step  was  to  find  the  value  of  y^,  the  impaction  efficiency, 

which  was  the  upper  limit  of  the  integral.  This  was  done  by  an 

iterative  procedure  called  the  half  interval  method  [16].  The  value  of 
y  for  the  critical  particle  was  estimated  far  upstream  and  the  path 
followed  to  the  target.  The  difference  between  the  ordinate  of  the 
tangent  path  and  the  ordinate  of  the  point  on  the  actual  path  parallel 
to  the  tangent  path  was  the  miss  criterion  used.  The  direction  of  the 

tangent  path  was  not  known  a  priori  buc  was  assumed  parallel  to  the 

actual  path.  The  half  interval  method  previously  mentioned  was  applied 
to  determine  a  better  initial  estimate.  Then  the  path  was  followed  to 
the  target  again  for  another  calculation  of  miss  distance.  This 
process  was  repeated  several  times  until  convergence  was  achieved.  The 
plane  of  initial  position  of  particles,  which  was  perpendicular  to  the 
flow  direction,  was  located  far  enough  from  the  target  so  that  free 
stream  conditions  prevailed.  A  distance  of  20  target  radii  upstream  of 
the  target  centre  was  found  to  be  sufficient.  The  particles  were 
assumed  to  travel  with  the  speed  and  direction  of  the  fluid  at  this 
upstream  position.  Starting  calculations  further  upstream  did  not 
change  results  significantly. 
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The  path  of  an  individual  particle  was  determined  step  by  step 

by  applying  a  fourth  order  Runge-Kutta  method  [16]  to  the  equations  of 

motion  (1)  and  (2).  The  value  K  in  these  equations  was  determined 

initially,  once  and  for  all,  for  the  whole  calculation  from  equation 

(4)  since  it  did  not  change  during  the  procedure.  The  value  of  Re  was 

calculated  for  each  time  step  using  equation  (6)  with  values  from  the 

previous  time  step  and  equation  (5)  in  which  ReQ  was  constant.  The 

value  of  CpRe  was  calculated  in  the  same  manner  by  numerical  solution 

of  equations  (11)  and  (12)  using  Newton's  method  [16]  for  finding  the 

zeros  of  a  function.  The  values  of  u  and  u  were  calculated  usinq 

x  y  y 

equations  (14a,  b). 

Once  the  value  of  y^  had  been  determined,  the  integration 
procedure  of  equation  (22)  was  applied  and  the  drag  coefficient  was 
calculated  from  equation  (20).  The  values  of  f(y)  used  in  equation 
(22)  were  determined  from  equation  (20)  for  each  particle  path  by 
tracing  a  particle  from  its  initial  position  to  the  cylindrical  target 
and  using  the  velocities  just  before  impaction  in  equations  (16)  and 
(17)  along  with  the  ordinate  of  the  impaction  point  in  equation  (18). 

Computer  programs  were  written  in  Fortran  for  the  DRES  Honeywell 
DP8/70C  using  the  CP-6  operating  system  to  calculate  impaction 
efficiency  and  drag  coefficient  for  various  particle  and  cylinder 
diameters.  The  complete  program  and  one  set  of  results  is  shown  in 
Annex  A.  The  application  of  the  method  is  straightforward,  with  the 
exception  of  the  step  by  step  integration  procedure  near  the  target 
cylinder.  This  procedure  is  therefore  explained  as  follows.  The 
particle  was  calculated  in  time  steps,  AT,  until  the  particle  was  one 
time  step  or  less  away  from  the  target.  Then  the  time  increment  size 
was  decreased  by  a  factor  of  ten  and  step  by  step  integration  continued 
until  the  particle  was  one  new  time  step  or  less  away  from  the  target. 
Finally,  the  time  increment  size  was  decreased  by  another  factor  of  ten 
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and  the  integration  allowed  to  proceed  to  the  target  again.  This 
method  ensured  that  the  position  of  the  target  and  particle  coincided 
within  a  maximum  error  given  by  the  distance  travelled  during  one 
percent  of  the  original  time  step  size,  while  allowing  the  particle  to 
reach  the  proximity  of  the  target  in  an  adequate  but  small  number  of 
steps,  to  use  computing  time  efficiently. 

RESULTS 


The  same  particle  tracing  technique  was  used  in  obtaining  the 
collection  efficiency  and  drag  coefficient.  Therefore,  a  comparison  of 
the  calculated  results  to  existing  information  on  collection  effici¬ 
ency,  for  which  previous  information  is  available,  is  helpful  to  gain 
confidence  in  the  results.  Collection  efficiencies  were  calculated  for 
a  wide  range  of  inertia  parameters  with  particle  specific  gravity  one 
(Table  1)  and  plotted  in  Figure  3  along  with  results  obtained  by  a 
mechanical  analogue  [14],  which  was  fundaments  1 ly  a  differential  ana¬ 
lyzer  constructed  for  solving  the  equations  of  motion  [17],  The  col¬ 
lection  efficiencies  were  compared  for  a  free  stream  Mach  number  of 
0.4  and  a  value  of  <t>  of  1000,  where  0  is  a  dimensionless  parameter  in¬ 
dependent  of  particle  size  given  by  the  following  equation: 

,  Re  2 
0  =  — o 

K  (23) 

The  air  properties  used  were  obtained  from  the  ICAO  Standard  Atmosphere 
at  an  altitude  of  3048  m  (10000  ft)  [18,  19].  The  compressible-air 
velocity  components  in  the  flow  field  from  equation  (14a,  b)  were  cal¬ 
culated  at  a  Mach  number  of  0.4  because  the  effect  of  compressibility 
for  completely  subsonic  flow  is  considered  greatest  at  or  near  the  free 
stream  critical  Mach  number  [14].  The  value  obtained  for  the  free 
stream  critical  Mach  number  depends  on  the  number  of  correction  terms 
applied  to  the  incompressible  flow  field  and  appears  to  approach  0.4  as 
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a  limit  [15].  The  value  of  1000  for  0  was  chosen  as  a  representative 
middle  value  from  the  results  reported  in  reference  [14]. 

Collection  efficiencies  were  also  calculated  for  incompressible 
flow  velocity  using  the  same  values  of  K,  0  and  atmospheric  conditions 
as  for  the  compressible  flow.  This  was  accomplished  simply  by  setting 
M  =  0  in  equations  (14a, b)  and  performing  the  calculations.  The  col¬ 
lection  efficiencies,  along  with  the  drag  coefficients  for  both  com¬ 
pressible  flow  at  Mach  number  0.4  and  incompressible  flow  are  shown  in 
Table  1.  The  collection  efficiencies  for  i ncompress ibl e  flow,  are  also 
plotted  in  Figure  4,  along  with  results  obtained  from  the  mechanical 
analogue  previously  described  [17,  21]. 

Calculation  of  collection  efficiencies  and  drag  coefficients  for 
a  wide  range  of  inertia  parameters,  using  particle  specific  gravity  of 
2.6,  the  density  of  sand,  and  dimensionless  parameters,  0,  were 
performed  using  the  ICAO  standard  atmosphere  at  sea  level  (Table  2). 
These  are  intended  mainly  to  provide  information  to  engineers  who  are 
interested  in  blast  hardening  of  structures,  but  they  could  also  be 
applied  to  estimate  collection  efficiencies  of  drops  or  solid 
particulates  in  chemical  defence  problems.  The  calculated  results  of 
collection  efficiencies  are  shown  along  with  fitted  curves  for  Mach 
number  0.4,  and  incompressible  flow,  in  Figures  5  and  6,  respectively. 
Similarily,  the  drag  coefficients  are  shown  in  Figures  7  and  8. 

An  example  of  the  distribution  of  forces  for  elastic  reflection 
of  particles  from  a  cylinder  is  shown  in  Figure  9,  where  the  function 
vxo^vx  -  v  cos  1  ^  Pitted  as  a  function  of  the  perpendicular 
distance  from  a  line  through  the  center  of  the  cylinder  along  the  free 
stream  flow  direction.  Comparison  of  surface  velocity  on  a  cylinder 
between  an  incompressible  flow  field  and  a  compressible-air  flow  field 
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at  Mach  number  0.4  is  shown  in  Figure  10  [14], 

DISCUSSION 


The  calculated  collection  efficiencies  agree  favourably  with  the 
theoretical  results  of  other  workers  for  which  a  mechanical  analogue 
was  used,  as  can  be  seen  in  Figures  3  and  4.  The  difference  in  results 
may  be  at  least  partly  explained  by  the  fact  that  the  equations  of 
motion,  (1)  and  (2),  were  considered  mathematically  in  two  different 
ways.  The  other  workers  [14,  17]  used  a  linearized  approximation  for 
the  equation  of  motion  between  x  =  -»  and  x  =  -5.  The  equations  were 
then  integrated  directly  from  x  =  -5  to  the  cylinder  using  the  starting 
conditions  at  x  =  -0.5  obtained  from  the  linearized  equations  [17], 
The  present  method  integrates  the  equations  of  motion  directly  from  x  = 
-20.0  to  the  cylinder.  The  results  agree  very  closely  for  the  highest 
values  of  the  inertia  parameter,  but  differ  more  for  the  smaller 
values.  Experimenting  with  various  time  step  sizes  showed  that  the 
results  were  sensitive  to  variation  in  this  parameter.  All  the  results 
shown  in  Figures  3  to  8  were  calculated  with  a  time  step  size  of  AT  = 
0.1,  which  was  chosen  such  that  convergence  was  provided  while  keeping 
computation  time  reasonably  short.  Using  this  time  step  size,  one  data 
point  could  be  obtained  in  30  minutes  from  an  IDM  RESEARCH  T286 
operating  at  11.7  mega-Herz,  with  math  co-proces sor ,  while  6  minutes 
were  required  with  the  ORES  Honeywell  DP8/70C  and  1  minute  with  the 
ORES  FPS  array  processor.  The  effect  of  reducing  AT  by  a  factor  of  5 
is  shown  for  one  value  of  #  and  seven  values  of  K  in  Table  3.  The  dif¬ 
ferences  due  to  the  time  step  difference  in  the  collection  efficiencies 
and  drag  coefficients  over  the  whole  range  of  inertia  parameters  are 
less  than  0.026  and  0.061,  respect i vely .  The  percentage  errors  in¬ 
crease  with  decreasing  values  of  K  due  to  the  fact  that  the  absolute 
values  decrease  with  K,  as  can  also  be  seen  in  Table  3.  The  effect  of 
varying  time  step  size  over  a  wide  range  is  shown  in  Table  4  for  the 
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smallest  value  of  K  considered.  The  results  indicate  that  the  value  of 
the  calculated  point  is  very  nearly  a  linear  function  of  time  step  size 
for  both  collection  efficiency  and  drag  coefficient.  The  intercept  for 
AT  =  0  is  0.0634  for  collection  efficiency,  which  falls  slightly  below 
the  results  presented  in  reference  [14],  as  can  be  seen  by  comparison 
to  the  curve  plotted  in  Figure  3.  The  fluid  velocity  components  given 
by  the  long  equations  (14a,  b)  were  plotted  as  shown  in  Figure  10.  The 
calculated  results  shown  agree  with  those  in  references  [14]  and  [15], 
which  contain  identical  figures,  as  a  comparison  check  showed. 

There  is  some  experimental  data  available  for  comparison  to 
calculated  theoretical  results.  Measurements  of  the  collection  effici¬ 
ency  of  cylinders  by  Ranz  and  Wong  [22]  represent  the  main  body  of  data 
on  which  comparisons  of  theory  and  experiment  are  based  [20].  The  ex¬ 
periments  were  conducted  with  a  3  mil  (76.2  pm)  platinum  wire  as  the 
cylindrical  collector  and  air  velocities  ranging  from  428  to  5090  cm 
sec'1.  Monodisperse  aerosols  of  sulfuric  acid  with  diameters  ranging 
from  0.56  to  1.30  pm  were  set  up  using  a  condensation  aerosol  genera¬ 
tor.  The  results  of  these  experiments  are  shown  in  Figure  10  [20], 
along  with  the  theoretical  curve  [21]  produced  by  the  methods  described 
in  reference  [17].  The  conditions  of  the  experiments  were  also  simula¬ 
ted  for  one  data  point  with  the  present  calculation  method,  and  the 
result  plotted  in  Figure  10.  This  calculation  applied  the  Cunningham 
correction  factor  for  small  particles  given  by  equation  (13).  The 
factor  for  the  selected  particle  size  of  1.073  pm  was  applied  to  the 
particle  drag  coefficient  in  air  for  all  calculations  to  obtain  the  re¬ 
sult.  The  Cunningham  correction  factor  was  not  applied  to  the  calcula¬ 
ted  results  shown  in  Figures  3  to  8,  but  the  effect  of  neglecting  it 
can  be  seen  by  examining  Table  5.  The  effect  of  not  applying  the  cor¬ 
rection  factor,  or,  in  other  words,  considering  the  fluid  of  continuum, 
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is  less  than  2  percent  for  both  the  collection  efficiency  and  drag  co¬ 
efficient  for  all  but  the  smallest  particle  sizes  considered.  This 
certainly  applies  to  the  portions  of  the  curves  represented  by  solid 
lines  for  compressible  flow  and  also  for  incompressible  flow  since  the 
same  size  particles  were  considered. 

As  mentioned  previously,  there  is  also  some  experimental  infor¬ 
mation  on  drag  forces  due  to  particles  in  air  flow  around  cylinders  [6, 
7],  These  experiments  were  performed  with  spherical  glass  beads  of 
0.0055,  0.0155  and  0.0470  pm  diameters  moving  vertically  toward  a 
cylinder  0.2  cm  diameter  in  a  tube.  The  drag  coefficients  for  the  par¬ 
ticles  were  adjusted  to  account  for  the  effect  of  gravity  in  retarding 
the  velocity  of  the  particles  [8].  The  drag  coefficients  obtained  were 
between  2.3  and  2.4  for  $  between  1.5  and  15  and  inertia  parameters 
ranging  from  18  to  1300,  which  is  to  be  compared  approximately  to  the 
upper  curve  in  the  high  range  of  inertia  parameters  shown  in  Figure  8. 

The  results  shown  in  Table  2  and  plotted  in  Figures  5  to  8  in¬ 
dicate  that  the  effect  of  compressibility  on  the  drag  due  to  particles 
large  enough  so  that  the  air  can  be  considered  a  continuum  is  negligi¬ 
ble  in  comparison  to  the  effect  of  compressibility  on  drag  due  to  the 
air  itself.  This  can  be  seen  by  comparison  to  the  results  shown  in 
references  [1,  2,  3,  4  and  5],  which  deal  with  drag  in  compressible 
flow.  The  reason  that  the  effect  of  compressibility  on  particle  tra¬ 
jectories  and  hence  collection  efficiency  and  drag  coefficient  is  neg¬ 
ligible  is  found  in  the  comparison  of  the  compressible  air  velocity 
field  with  the  incompresible  air  velocity  field.  The  effect  of  compres¬ 
sibility  is  most  pronounced  close  to  the  cylinder  [14].  Furthermore, 
the  effect  of  compressibility  is  pronounced  only  over  a  small  portion 
of  the  cylinder  near  the  top  and  bottom,  as  shown  in  Figure  10. 
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CONCLUSIONS 


The  collection  efficiences  and  drag  coefficients  of  cylinders 
due  to  airborne  particles  have  been  calculated  over  a  wide  range  of 
inertia  parameters  and  Reynolds  numbers  for  incompressible  flow  and 
subsonic  compressible  flow  near  critical  Mach  number.  The  results 
agree  favourably  with  the  small  amount  of  experimental  data  available. 
The  effect  of  compressibility  on  the  drag  forces  due  to  particles 
carried  in  an  air  flow  was  found  negligible  for  particles  large  enough 
so  that  the  fluid  can  be  considered  a  continuum. 
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TABLE  1 

CALCULATED  COLLECTION  EFFICIENCIES  AND  DRAG  COEFFICIENTS 
FOR  MACH  NUMBER  0.4  IN  STANDARD  ICAO  ATMOSPHERE 
AT  3048  METRES  HEIGHT 


INERTIA 

PARAMETER 

K 

CYLINDER 

DIAMETER 

D  cm 

PARTICLE 

DIAMETER 

d  pm 

DIMENSIONLESS 

PARAMETER 

0 

M  = 

=  0.4 

INCOMPRESSIBLE 

E 

m 

CD 

Em 

CD 

0.4 

1.631 

1000 

0.0913 

0.0357 

0.1040 

0.0490 

1 

1.631 

1000 

0.2621 

0.2670 

0.2794 

0.3118 

4 

1.631 

1000 

0.5728 

1.0798 

0.5896 

1.1483 

10 

1.631 

1.356 

1000 

0.7404 

1.6536 

0.7536 

1.7124 

40 

1.631 

2.730 

1000 

0.8924 

2.2292 

0.8990 

100 

1.631 

4.317 

1000 

0.9410 

2.4212 

0.9449 

2.4401 

400 

1.631 

8.634 

1000 

0.9747 

2.5558 

0.9763 

2.5649 

UNCLASSIFIED 


UNCLASSIFIED 


TABLE  2 


CALCULATED  COLLECTION  EFFICIENCIES  AND  DRAG  COEFFICIENTS 
FOR  COMPRESIBLE  FLOW  AT  MACH  NUMBER  0.4  AND 
INCOMPRESSIBLE  FLOW  IN  STANDARD  ICAO  ATMOSPHERE  AT  SEA  LEVEL 


INERTIA 

CYLINDER 

PARTICLE 

DIMENSIONLESS 

M  = 

0.4 

INCOMPRESSIBLE 

PARAMETER 

K 

DIAMETER 

D  cm 

DIAMETER 
d  urn 

PARAMETER 

E 

C 

E 

C 

m 

D 

m 

D 

0.4 

0.02531 

0.215 

10 

0.1384 

0.0811 

0.1533 

0.1042 

1 

0.02531 

0.339 

10 

0.3644 

0.5040 

0.3807 

0.5588 

4 

0.02531 

0.679 

10 

0.6999 

1.5311 

0.7110 

1.5823 

10 

0.02531 

1.073 

10 

0.8422 

2.0471 

0.8489 

2.0791 

40 

0.02531 

2.146 

10 

0.9454 

2.4416 

0.9479 

2.4554 

100 

0.02531 

3.393 

10 

0.9721 

2.5465 

0.9733 

2.5545 

400 

0.02531 

6.787 

10 

0.9878 

2.6092 

0.9882 

2.6133 

0.4 

0.2531 

0.679 

100 

0.1190 

0.0611 

0.1330 

0.0802 

1 

0.2531 

1.073 

100 

0.3283 

0.4130 

0.3451 

0.4665 

4 

0.2531 

2.146 

100 

0.6532 

1.3589 

0.6665 

1.4179 

10 

0.2531 

3.393 

100 

0.8056 

1.9037 

0.8151 

1.9480 

40 

0.2531 

6.787 

100 

0.9272 

2.3691 

0.9312 

2 . 3889 

100 

0.2531 

10.731 

100 

0.9618 

2.5051 

0.9639 

2.5168 

400 

0.2531 

21.462 

100 

0.9837 

2.5924 

0.9845 

2.5981 

0.4 

2.531 

2.146 

10UU 

0.0913 

0.0357 

0.1040 

0.0490 

1 

2.531 

3.393 

1000 

0.2621 

0.2671 

0.2794 

0.3118 

4 

2.531 

6.787 

1000 

0.5729 

1.0772 

0.5897 

1.1485 

10 

2.531 

10.731 

1000 

0.7405 

1.6540 

0.7536 

1.7127 

40 

2.531 

21.462 

1000 

0.8925 

2.2294 

0.8991 

2.2609 

100 

2.531 

33.93 

1000 

0.9410 

2.4213 

0.9449 

2.4401 

400 

2.531 

67.87 

1000 

0.9747 

2.5558 

0.9763 

2.5649 

0.4 

25.31 

6.787 

10000 

0.0525 

0.0122 

0.0622 

0.0183 

1 

25.31 

10.731 

10000 

0.1779 

0.1260 

0.1943 

0.1567 

4 

25.31 

21.462 

10000 

0.4556 

0.7100 

0.4747 

0.7823 

10 

25.31 

33.934 

10000 

0.6341 

1.2707 

0.6513 

1.3401 

40 

25.31 

67.867 

10000 

0.8249 

1.9596 

0.8362 

2. 0107 

100 

25.31 

10/. 308 

1000U 

0.8970 

2.2426 

0.9043 

2.2764 

400 

25.31 

214.615 

10000 

0.9531 

2.4671 

0.9566 

2.4847 

0.4 

126.57 

15.176 

50000 

0.0306 

0.0044 

0.0378 

0.0072 

1 

126.57 

23.995 

50000 

0.1191 

0.0592 

0.1319 

0.0759 

4 

126.57 

47.989 

50000 

0.3537 

0.4464 

0.3733 

0.5065 

10 

126.57 

75.876 

50000 

0.5284 

0.9099 

0.5484 

0.9892 

40 

126.57 

151.756 

50000 

0.7441 

1.6431 

0.7594 

1.7125 

100 

126.57 

239.947 

50000 

0.8340 

2.0057 

0.8494 

2.0551 

400 

126.57 

479.894 

50000 

0.9197 

2.3288 

0.9262 

2.3592 

UNCLASSIFIED 


UNCLASSIFIED 

TABLE  3 

EFFECT  OF  TIME  STEP  INCREMENT  SIZE 
ON  COLLECTION  EFFICIENCY  AND 
DRAG  COEFFICIENT  FOR  VARIOUS  INERTIA  PARAMETERS 


ICAO  ATMOSPHERE  AT  3048  m  M=0.4  K=0.4  *  =  1000 


K 

T) 

RATIO 

CD 

RATIO 

at=o.i 

AT=0.02 

m  AT=0.1 

AT=0 . 1 

AT=0.02 

cdat=o.i 

E 

m  AT=0.2 

CpAT-0.02 

0.4 

0.09127 

0.06882 

1.3262 

0.0357 

0.0210 

1.7000 

1 

0.2621 

0.2361 

1.1101 

0.2670 

0.2194 

1.2170 

4 

0.6728 

0.5553 

1.0315 

1.0768 

1.0158 

1.0601 

10 

0.7404 

0.7296 

1.0148 

1.6536 

1.6077 

1.0286 

40 

0.8924 

0.8879 

1.0051 

2.2292 

2.2069 

1.0101 

100 

0.9410 

0.9386 

1.0026 

2.4212 

2.4087 

1.0052 

400 

0.9742 

0.9737 

1.0005 

2.5558 

2.5507 

1.0020 

UNCLASSIFIED 


UNCLASSIFIED 


TABLE  4 

EFFECT  OF  TIME  STEP  INCREMENT  SIZE  ON 
CALCULATED  VALUES  OF  COLLECTION  EFFICIENCY  AND 
DRAG  COEFFICIENT  FOR  SMALL  INERTIA  PARAMETER 


ICAO  ATMOSPHERE  AT  3048  m  M  =  0.4  K  =  0.4  *=  1000 


AT 

COLLECTION  EFFICIENCY 

DRAG  COEFFICIENT 

CALCULATED 

POINT 

FITTED 

LINE 

CALCULATED 

POINT 

FITTED 

LINE 

0.1 

0.0913 

0.0912 

0.0357 

0.0354 

0.05 

0.0771 

0.0773 

0.0260 

0.0265 

0.02 

0.0688 

0.0690 

0.0210 

0.0212 

0.004 

0.064G 

0.0645 

0.0186 

0.0184 

0.001 

0.0638 

0.0637 

0.0181 

0.0179 

FITTED  LINE 

r2 

0.99981 

0.99750 

y  =  a  +  bx 

a 

0.006340 

0.01768 

r2  =  coefficient  of 

b 

0.2779 

0.1772 

determination 

UNCLASSIFIED 


UNCLASSIFIED 

TABLE  5 

EFFECT  OF  CUNNINGHAM  CORRECTION  FACTOR 
FOR  VERY  SMALL  PARTICLES  ON  COLLECTION  EFFICIENCY  AND 
DRAG  COEFFICIENT  AT  MACH  0.4 


PARTICLE 

CUNNINGHAM 

INERTIA 

DIMENSIONLESS 

CORRECTED 

DIAMETER 

CORRECTION 

PARAMETER 

PARAMETER  0 

UNCORRECTED 

d  pm 

FACTOR  C 

c 

K 

0 

Emc/E 

m 

CDc/ 

ld 

1.073 

1.153 

10 

10 

1.014 

m 

2.146 

1.076 

40 

10 

1.001 

HU 

2.146 

1.076 

4 

100 

1.011 

2.146 

1.076 

0.4 

1000 

1.106 

1.275 

3.393 

1.048 

1 

1000 

1.008 

1.010 

6.787 

1.024 

0.4 

10000 

1.002 

1.008 

UNCLASSIFIED 
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SM  1267 


UNCLASSIFIED 


CIRCULAR  CYLINDER 


VELOCITY  CHANGE  OF  A  PARTICLE  DUE  TO  INTERACTION  WITH  A  CYLINDER 


UNCLASSIFIED 


CYLINDER  COLLECTION  EFFICIENCY  IN  COMPRESSIBLE  AIR  FLOW  AT  MACH  NUMBER 
OF  0.4  AND  0  OF  1000  IN  ICAO  STANDARD  ATMOSPHERE  AT  3048  METERS 


UNCLASSIFIED 


IIM  ICAO  STANDARD  ATMOSPHERE  AT  3048  METERS. 


UNCLASSIFIED 


SM  1267 


Figure  7 

CALCULATED  DRAG  COEFFICIENT  FOR  CYLINDER  DUE  TO  PARTICLES  ALONE  AT 
MACH  NUMBER  OF  0.4  IN  ICAO  STANDARD  ATMOSPHERE  AT  SEA  LEVEL 


UNCLASSIFIED 


DISTRIBUTION  FUNCTION  OF  DRAG  COEFFICIENT 


UNCLASSIFIED 


SM  1267 


DISTANCE  FROM  SYMMETRICAL  CENTRE  LINE  Y 


Figure  9 

DISTRIBUTION  OF  FORCE  DUE  TO  ELASTIC  REFLECTION  OF  PARTICLES  FROM  A 
CYLINDER  AT  MACH  NO.  OF  0.4  in  ICAO  STANDARD  ATMOSHERE  AT  SEA  LEVEL 


UNCLASSIFIED 


TANGENTIAL  SURFACE  VELOCITY,  us,  RATIO  TO  FREE  STREAM  VELOCITY 


UNCLASSIFIED 


SM  1267 


CYLINDER  ANGLE,  d 


Figure  10 

COMPARISON  OF  SURFACE  VELOCITY  ON  A  CYLINDER  BETWEEN 
INCOMPRESSIBLE  FLOW  FIELD  AND  COMPRESSIBLE-AIR  FLOW  FIELD  AT  MACH  NO.  0.4 


UNCLASSIFIED 


COLLECTION  EFFICIENCY  E 


UNCLASSIFIED 


SM  1267 


Figure  11 

COMPARISON  OF  PRESENT  CALCULATION  METHOD  TO  THEORY  AND 
EXPERIMENT  OF  OTHER  INVESTIGATORS 
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